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We present a solution of Kitaev's spin model on the honeycomb lattice and of related topologically ordered 
spin models. We employ a Jordan- Wigner type fermionization and find that the Hamiltonian takes a BCS type 
form, allowing the system to be solved by Bogoliubov transformation. Our fermionization does not employ 
non-physical auxiliary degrees of freedom and the eigenstates we obtain are completely explicit in terms of the 
spin variables. The ground-state is obtained as a BCS condensate of fermion pairs over a vacuum state which 
corresponds to the toric code state with the same vorticity. We show in detail how to calculate all eigenstates 
and eigenvalues of the model on the torus. In particular, we find that the topological degeneracy on the torus 
descends directly from that of the toric code, which now supplies four vacua for the fermions, one for each 
choice of periodic vs. anti-periodic boundary conditions. The reduction of the degeneracy in the non-Abelian 
phase of the model is seen to be due to the vanishing of one of the corresponding candidate BCS ground-states 
in that phase. This occurs in particular in the fully periodic vortex-free sector. The true ground-state in this 
sector is exhibited and shown to be gapped away from the three partially anti-periodic ground-states whenever 
the non-Abelian phase is gapped. 
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I. INTRODUCTION 

A combination of special properties has made Kitaev's spin 
model on the honeycomb lattice^ a very popular subject of 
study in recent years. The model has a basic Hamiltonian 
with only nearest neighbour interactions, but nevertheless, ex- 
tensions of the model with magnetic field like terms have 
both Abelian and non-Abelian topological phases. Moreover 
these extended models can still be exactly solvable, allowing 
in principle for direct study of both phases and of the phase 
transition. 

In spite of the availability of what is by now a large col- 
lection of exact solutiona 'i 2 ! 3 ! 4 ! 5 ! 6 ! 7 ! 8 ! 9 ! 10 ] 1 'i 12 i 13 i 14 i 15 , all based 
on fermionization techniques, the two types of phases of the 
system are generally understood through the use of different 
methods of analysis. On the one hand Kitaev showed 1 using 
perturbation theory that the low energy effective theory of the 
Abelian phases is equivalent to his Z2 toric code modeli£ and 
in the mean time, extensive further perturbative work on this 
phase has been don e 17 i 18 i 19 i 20 i 21 . On the other hand, the non- 
Abelian phase is only understood using the fermionized exact 
solutions of the system. Its topological order is known to be 
described by the Ising model of topological field theory. 

There have been a number of recent works linking the Z2 
toric code and Ising topological field theories directly. For ex- 
ample it was demonstrated that a Z2 toric code theory could be 
formed by condensing bosonic excitations in a doubled Ising 
theory 22,23 . Also, excitations with properties of Ising anyons 
were constructed from superpositions of the electric and mag- 
netic excitations of the toric code model^ 4 .. In the present 
work, we want to explore the relationship between these two 
types of topological order within the context of the honey- 
comb model. 

In order to do this, we introduce yet another solution of 
the model, but one which is particularly useful for studying 
the relation between the toric code and Ising type topologi- 



cal orders which exist in the model. Our solution is again 
by a Jordan-Wigner type fermionization, however, the Jordan- 
Wigner transformation we employ is closely linked to a choice 
of basis for the Hilbert space adapted to perturbative analysis 
of the Abelian, toric code type phas o 17 ' 18 ' 19 . The fermions 
we use are also closely related to the deconfined fermionic 
excitations which were shown to occur throughout the phase 
diagram in ref.H^and which correspond to the fermionic exci- 
tations of the toric code in the Abelian phase. After fermion- 
ization, the model can be solved exactly and, as with other 
fermionization methods (cf. refs. HUD), the ground-state sec- 
tor of the system can be transformed to that of a spinless 
p-wave superconductor, as analyzed by Read and Green in 
ref.|2(j The ground-state is thus a BCS type state 27 , and can 
be related to the v = 5/2 fractional quantum Hall state of 
Moore and Read 28 . With our method, we obtain a vacuum 
for the fermionized theory which is exactly defined in terms 
of toric code stabilizers and independent of the couplings of 
the model. The ground-state for the full system, valid for all 
parameter space, is in fact a BCS type condensate over the 
toric code ground-state. Because the vacuum is independent 
of the coupling parameters, the mechanism for switching be- 
tween topological phases is contained exclusively within the 
BCS product. On the other hand, the topological degeneracies 
of the model are already present at the level of the toric code 
vacuum. The BCS product only lifts some of this degeneracy 
in the non-Abelian phase. 

The structure of the paper is as follows: We start with a 
short review of the model and of the emergence of the toric 
code as an effective description of the Abelian phase. We give 
special attention to a description of the Hilbert space of the 
model in terms of hard core bosons and effective spins on a 
square lattice, as this is essential preparation for our fermion- 
ization scheme. In section [Til] we fermionize and solve the 
model on the plane and give an explicit expression for the 
ground-state of the model in eq. d27l i. This expression in- 
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volves only the physical degrees of freedom of the model; no 
auxiliary variables are introduced anywhere in this work. In 
section llVl we extend our fermionization method to the torus. 
We explain how to construct the eigenstates in any given vor- 
ticity sector and how to calculate their corresponding energy 
eigenvalues. Since the creation of an odd number of fermions 
does not preserve vorticity, certain low energy eigenvalues 
which might be expected do not occur in every vorticity sec- 
tor. We give particular attention to the vortex free sector which 
contains the model's ground-states, and show that the energy 
of the lowest lying fully periodic state in this sector is lifted 
in the non-Abelian phase, proving that the ground-state of the 
non-Abelian phase is three-fold degenerate. In appendix [A] 
we give a general discussion of Hartree-Fock-Bogoliubov the- 
ory with gauge violating fermions, which is used as back- 
ground for the discussion in sectionllVl 

II. HONEYCOMB MODEL AND TORIC CODE 

A. Spin Hamiltonian and loop symmetries 

The system consists of spins on the sites of a hexagonal 
lattice. The Hamiltonian can be written as 

H = - (!) 

aE{x,y,z} i,j 

where Kg- = (rfcr? denotes a directional spin exchange in- 
teraction occurring between the sites i,j connected by a a- 
link see FIG. Q] We define a the basic unit cell of the lat- 
tice with the two unit vectors n x and n y as shown in FIG. 
[TJ By contracting each z-link to a single point we define 
the position vector labeling the z-dimers on a square lattice 
as q = q x n x + q y n y . 

Consider now products of K operators along loops on the 
lattice , Kf!-™ Kf™ ....K% ln) , where a (m ^ € x,y,z. Any 
loop constructed in this way commutes with the Hamiltonian 
and with all other loops. The shortest such loop symmetries 
are the plaquette operators 

W q = ululalalutul, (2) 

where the numbers 1 through 6 label lattice sites on single 
hexagonal plaquette, see FIG. Q] We will use the convention 
that q denotes the z-dimer directly below the plaquette. The 
fact that the Hamiltonian commutes with all plaquette opera- 
tors implies that we may choose energy eigenvectors | n) such 
that W q = (n \W q \ n) = ±1. If W q = -1 then we say that 
the state | n) carries a vortex at q. When we refer to a partic- 
ular vortex-sector we mean the subspace of the system with a 
particular configuration of vortices. The vortex-free sector for 
example is the subspace spanned by all eigenvectors such that 
W q = 1 for all q. 

On a torus, the plaquette operators are not independent, as 
they obey Y[ Wq = I where the product is over all q. There 
are also two independent homologically non-trivial loop sym- 
metries. To represent these we are free to choose any two 
closed loop operators that traverse the torus as long as they 
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FIG. 1: The plaquette operator W and the fermionic string S 

cannot be deformed into each other by plaquette multiplica- 
tion. All other homologically non-trivial loop symmetries can 
be constructed from the products of these two operators and 
the N/2 — 1 independent plaquette operators, (cf. ref. l20h . 
Note that when the torus is specified by periodic boundary 
vectors (x, y) which are integer multiples of the unit vectors 

i.e. x = N x n x and y = N y n y , it is natural to use overlapping 

i x \ i i 

products of alternating z- and x-links (Lq y = WK^K^) 

and alternating z- and y-links (Lq V J = Y[ K^K v ^ k ), as homo- 
logically non-trivial symmetries. We will generally use the 

(x) (v) 

operators L and L that run through the origin as the two 
independent symmetries. 

The model contains 4 distinct phases 1 . There are three 
gapped phases (1) A x with J x > J y + J z , (2) A y with 
J y > Jx + Jz and (3) A z with J z > J x + J y with (4) a 
gapless B phase existing in the parameter space between the 
three A phases. As each A phase is related to the others by a 
lattice rotation we confine our analysis to the A z and B phases 
with out loss of generality. 

The Hamiltonian (Q~|i is often extended to include perturbing 
terms Hi that (i) are sums of K operator products (ii) open 
a gap in the B-phase (iii) break time-reversal symmetry (T- 
symmetry), se o 1 ' 2 ' 3 ' 4,5 and the general analysis of the link or 
bond algebras in 14 . The breaking of T-symmetry is essential 
for relating the model to chiral p-wave superconductors. As 
the procedure we will outline here gives the same physical 
results as the quoted references for generalised T-symmetry 
breaking we restrict the explicit calculations to the three-body 
term studied in references[l||2j|3|: 

q 1=1 

with the second summation running over the six terms 

]T P(q)W = of4<4 + ala\a% + (4) 
i=i 

a\alal + afa^a z 6 + (rfafaf + ^of of . 

Recently a third type of phase has been discovered 
in extended honeycomb models, featuring gapped Dirac 
fermions 25 . The phase is opened by allowing the J a couplings 
to vary periodically on the lattice. Although we do not exam- 
ine this phase in this paper, the methodology employed below 
can in principle be used. 
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B. The Toric Code as an effective system 



C. Hard-core bosons and stabilizers 



The Hamiltonian ([T]) can be written in terms of hard- 
core bosons and effective spins of the z-dimers using the 
mapping-^: 



I T.TJ = I to), 
T.i D ) = I it A), 



I !■!□) = I 4,0), 
I !■!□) = ! 



(5) 



The labels on the left hand side indicate the states of the z- 
dimer in the S z basis. The first quantum number of the kets on 
the right hand side represents the effective spin of the square 
lattice and the second is the bosonic occupation number. The 
presence of a boson indicates an anti-ferromagnetic configu- 
ration of the spins connected by a z-link. 

In the A z -phase, the dominance of the J z means that spins 
on a z-dimer tend to align in the same direction, and there- 
fore in this limit the presence of bosons is energetically sup- 
pressed. A perturbative analysis for the low energy effective 
Hamiltonian in this regime shows that the first non-constant 
term, occurring at the 4th order, is 



HtC = -JeK^2Q q ®I 



(6) 



In the lattice orientation we use, r| changes the eigenvalues 
Qq-n x and Qq-n y while r| changes the eigenvalues Q q and 
Qq-n x -n y - On even-even lattices, which can be bi-colored, 
the quasi-particle excitations occur in two types usually la- 
beled e and m. Excitations of the same type are mutually 
bosonic but excitations of different types display Abelian any- 
onic statistics. Importantly, pairs of e and m particles behave 
as fermions. 

In the language of the stabilizer formalism , se o 37 ' 38 , 
we say that the TC states are stabilized by the operators 
Qq' I {Qq}) — Qq' I {Qq})- On a plane for example, the TC 
ground-state is the state such that Q q ,\ {Q q }) = \ {Qq}) for 
all q. However, despite their simple description in terms of 
the stabilizer formalism, it is important to recognize that the 
TC states are structurally non-trivial and display unusual en- 
tanglement and geometric properties. It has been shown for 
example that the ground-state of the TC system is a Projected 
Entangled Pair State (PEPS) with virtual dimension D=2 (see 
for example ref. 1401) . and can also be described in terms of 
string-net condensates and loop model s 41 ' 42 . 



The basis ([5]) also describes anti-ferromagnetic configura- 
tions of the z-dimers through the bosonic occupation number 
and forms an orthonormal basis for the full honeycomb sys- 
tem. The Pauli operators of the origina l spin Hamiltonian can 
be written as (see r eferences 17,18fll) : 



= r x q {b\ 



bq) 



M - b q ) 



r z q {I 



2b q b q ) 



(7) 



where b< and b are the creation and annihilation operators for 
the hard-core bosons. In this representation the Hamiltonian 
(Q~|i becomes 



Ho = -Jx^l + KK+nM 



+n x + b q+rix ) 



Jy^iT Z q {b\ - b q )Tl +ny {b\ + ny +b q+nv ) 
1 

j z J2(i-HK)- 



(8) 



with Q q = T z q Tl +n j v q+ny T* +n where is the Pauli op- 
erator acting on the effective spin at position q and J e ff = 

J 2 J 2 

ie 7j v i 3 — ■ This effective Hamiltonian, defined now on a square 
lattice, is unitarily equivalent to what is known as the toric 
code (TCA 

The operators Q q , like the plaquette operators W q , all 
commute with each other. The eigenvalues of each oper- 
ator can therefore be used as quantum numbers to specify 
eigenstates of the system. We write | {Q q }) where {Q q } is 
a full list of Q q eigenvalues. Excitations of the TC system 
are made/moved by applying r 2 and/or t v operators to a site. wnere jy 



with the perturbative term Hi given pictorially in FIG. [3] 
The Hamiltonian ([8]) has been used in the gapped A z phase 
to perturbatively calculate effective Hamiltonians and other 
measures, to the 10th order in some case a 17 ' 18 ' 19 . In the next 
section we will show how to fermionize this Hamiltonian by 
attaching string operators to the hard-core bosons. The pro- 
cedure is much like other Jordan- Wigner type approaches but 
the operator strings that we choose, will be tailored for this 
system. 

In this representation the plaquette operators (O become 



W q = (I-2N q )(I~2N q+riy )Q q 



(9) 



i ~ ^</^V Thi s relation (O is very useful be- 
cause it allows one to write down an orthonormal basis for 
the full honeycomb system^. The basis can be written as 
I {Wg}) {?}), where {q} lists the sites with non-zero bosonic 
occupation and the eigenvalues {M^q} determine the vortex 
sector. Note that to determine the structure of the state one 
still uses the effective operators Q q as the stabilizers with 
the eigenvalues reflecting the vorticity {Wq} through ([9). In 
the special case where there are no broken dimers we have 
| {Wq}, {0}) = | {Qq}) with W q = Qq for all q. 

To specify a state on a torus we also give two additional 
quantum numbers associated with the homologically non- 
trivial loop symmetries. We can normally choose these to be 

(x) (y) (x) 

the eigenvalues Zq and l of the independent operators L 
and Lq J ^ described above and the generic state in this case can 

be written as | {W q }, {<?}) where the list {W q } con- 

tains N x N y — 1 independent Wq's. 

In the next section we will have cause to use a generali- 
sation of the expression (O for products of plaquette oper- 
ators. Of particular importance, because of the conventions 
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used, will be the products arranged vertically on the effective 
lattice. We have in this case 



9y-l 

q' v =o 



y 



(10) 



1y-l 



= (I-2N q „, )(I-2N qxAy ) H Q q ^ q , 

l'y=° 

and we see that only the bosonic occupation numbers at the 
upper and lower left corners of the plaquette product need to 
be taken into account. 
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FIG. 2: Bosonic and effective spin decomposition of the operator 
string S. 



III. FERMIONIZATION 

We now show how to turn the hardcore bosons into 
fermions using a new Jordan- Wigner type transformation that 
is designed for this model and in particular for the basis 
I {Wg}i {?}) described above. The procedure has anumberof 
advantages over other fermionization techniques. For exam- 
ple, the method does not introduce additional un-physical de- 
grees of freedom like the Majorana approach originally used 
to solve the problem 1 but still allows a transparent encoding 
of the vorticity within the fermionic Hamiltonian. In addi- 
tion the procedure also reveals much more about the actual 
eigenstates of the system than previous fermionization meth- 
ods. We will see for example, like Chen and Nussinov-, that 
the ground-state of the system is a BCS type product acting on 
the vacuum. However, our vacuum will be a exactly defined in 
terms of Toric Code stabilizers meaning that the eigenstates of 
the system can be written in simple closed form expressions 
that do not require implicit spectral projection. 

We begin by defining a particular string operator using 
overlapping products of the Kf^ terms of the original Hamil- 
tonian. The string will serve two purposes: (i) it will break/fix 
z-dimers at a single location q of the lattice thereby creat- 
ing/annihilating hard-core bosons (ii) it will enforce inter-site 
fermionic commutation relations effectively turning our hard- 
core bosons into fermions. 

Our convention will be to first apply a single a x term to a 
black site of the z-link which we set to be the origin. The rest 
of the string is made by applying first alternating Kfj and KJ k 
until we reach a required length and then apply alternating 
Kf m and K^ n terms ending on the black site at q, see FIG.Q] 
Explicitly we write 



S a 



y y % 

CT (9*,?y),B%a ; ,? y -i),n fT (<fc,<Z!/-i)n 

y y z z x 

•'•%x,l),B fT (9 ;E ,0), n Cr (fe,0), D cr (g x ,0) B o '( 9:c ,0), 



(11) 
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l <7 (0,0) n cr (0,0),p°'(0,0) ; 



'(0,0), 



Using the representation o f 17 ' 18 ' 19 we can decompose dTTb into 
the effective spin and bosonic subspaces, i.e. S = S e <g> Sb- 
In this decomposition there are four different types of struc- 
tures to observe on the effective lattice: (1) the line including 
the starting point A up to, but not including, the turning point 



B, (2) the turning point B = (q Xl 0), (3) the exclusive inter- 
val BC, and (4) the end point C — (q x , q y ), see FIG [2] and 
TABLE U 
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TABLE I: The string S as four unique segments. While bosons are 
only created/destroyed at the endpoint C of the string, the sites in the 
[A, B) interval also have non-trivial bosonic dependence. 

The operator S q squares to unity while different operators 
S q , S q ' anti-commute with each other. This leads us to iden- 
tify the string S q with the following sum of fermionic creation 
and annihilation operators: S q = c q +c q = (b q +b q )S q where 

S q is simply the string S q but with the bosonic dependence 
of the end-point C removed, see TABLE U Individually our 
fermionic canonical creation and annihilation operators are 



C q — bqSq, C q — b q S q 



(12) 



where the strings now insure that the operators c q and c q obey 
the canonical fermionic anti-commutator relations 

H,c qt } = d gql , {ct,^,} = 0, {c q ,c gl } = 0. (13) 

The operators c q and c q must both create/annihilate vortices 
at — n x and — n x — n y . We can therefore think of each fermion 
as being bound to a vortex pair at the origin. Remarkably, this 
vortex-pair can itself be thought of as a fermion and can even 
be moved without changing the energy of the system 2 -. An- 
other interesting insight can be obtained by noting that the 
strings S' q , which we attach to the hard-core bosons to make 
them fermionic, change the eigenvalues of the Q q operators 
at — n x , —n x — n y , q and q — n y . The creation of the (A 
fermions is thus reflected in the effective spins by the creation 
of two bound fermionic e — to pairs. These associated states 
are structurally equivalent, up to the mapping ((5]), to the cor- 
responding TC states. 
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FIG. 3: The Kitaev three body term in the effective-spin/hard-core boson and fermionic notation. On a torus the values of X and Y are as 
those for the basic two-body Hamiltonian described below, except for the terms P(q)^ and P(q)^ when q y = N y . 



If we invert ( fl2l i and substitute the relevant expressions into 
the Hamiltonian dS) we get 



q+riy 



q+n.. 



-~q+n y 



.) 



+ Jz 



q 



(2c q c q 



I), 



(14) 



where, in the plane, Y q — I for all q and X q is defined in ( fTob - 
The fermionic representation of the perturbative term i?i is 
given in FIG. [3] We restrict the Hilbert space to the relevant 
vortex-configuration by replacing X q by the eigenvalues X q 
of that configuration. In the simplest case of the vortex free 
sector we have X q = 1 for all q. This sector, because of a 
theorem by Lieb 29 , is known to contain the system ground- 
state and can be solved exactly in the thermodynamic limit 
by moving to the momentum representation with the Fourier 
transform 



Jkq 



<v"-. (15) 
After substitution into ( fT4l i and anti-symmetrization we have 



tkc{c k + ~(A4.cL fe + A*c_ fc c fc ) 



1 



where 



with 



/'• 

0k 



A fc = a k + i(3 k 

-2J Z 

2J X cos(k x ) + 2J y cos(fc y ) 
AK{sm{k x ) — sm(ky) — sin(fc 2 
2J X sm(k x ) + 23 y sin(fc y ). 



-MJ Z (16) 



(17) 
(18) 

(19) 
(20) 

ky)) (21) 

(22) 



where the extended Hamiltonian Hi is now fully contained 
within in the ak term. The procedure also gives agreement 
with the other fermionization techniques to analyse the ex- 
tended mode l 2 ' 3 ' 4,7 ' 9 . We note in particular that the technique 
can be used to replicate the dispersion relations of 5 where the 
p-wave pairing can be tuned to have k x +ik y chiral symmetry 
thus allowing a direct link with the work of Read and Green 26 
and subsequent analysi s 30 ' 31 ' 32 ' 33 , relating the Pfaffian Quan- 
tum Hall states, p-wave superconductors and the Ising CFT 
model. 

The Hamiltonian ( fToT ) is diagonalized by Bogoliubov trans- 
formation 



Ik = u k c k - 

where u k and v k satisfy \u k \ 2 - 
H = EE k hhk- 1/2), with 



v k c 



k- 



(23) 



\v k \ 2 = 1. We then have 



E k 

Uk 



e k + \^\ 



= y/l/2(l + £ k /E k ) 
= iy/l/2(l-£ k /E k ) 



(24) 

(25) 
(26) 



The ground-state, annihilated by all 7^, and of energy E gs — 
~ 1 1 Ekdk, can be seen to be the BCS type state 

I gs) = ]J(u k + v k clcl k )\ {W q }, {0}). (27) 



This expression is, to the best of our knowledge, the first 
closed form expression of the ground-state that does not re- 
quire additional spectral projection. It is noteworthy because 
it combines two powerful wavefunction descriptors i.e. the 
BCS product and the Stabilizer formalism. In the expression, 
which is valid everywhere in the model's parameter space, 
the fermionic vacuum is fixed to be the toric code ground- 
state. While this implies that any mechanism for switching be- 
tween the Abelian and non-Abelian topological phases must 
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be contained exclusively within the BCS product, we should 
also recognize that the Abelian phase is Z 2 x Z 2 because the 
fermionic vacuum is Z 2 x Z 2 and not because of any myste- 
rious property of the BCS product. To see this more clearly 
note that in the A z phase with (J z = l,J x ,J y — > 0) we have 
Uk — » 1 and Vk — > and the the ground-state of the full sys- 
tem |gs) -» |{WJ,{0}) = \{Q q }), where Q q = 1 and 
W q = 1 for all q. This is of course what one expects from the 
perturbation theory, see for exampl&22i21. 



IV. FERMIONIZATION ON A TORUS 

The fermionization procedure above may be extended to 
systems that live on a torus. Going to the torus allows for 
the study of finite size systems without fixing boundary con- 
ditions. It also allows us to probe the topological order of 
the model's A and B-phases directly. From the predictions 
of topological-QFT the ground-state degeneracy on a torus 
should be equal to the number of topological sectors, or quasi- 
particle types in the system. In the Abelian A-phase we should 
have a 4-fold degeneracy because we have 4 distinct parti- 
cle types : the trivial particle or vacuum, the e particle, the 
m particle, and the fermionic e-m composite, which as we 
have mentioned, corresponds to our excitation. In the non- 
Abelian phase the theory predicts we have three distinct parti- 
cle types with the distinction between e and m particle types 
no longer applicable 1 . However, as far as we know, a direct 
analysis of the ground-state degeneracy in this phase has not 
been done. This may be because on the surface, it appears that 
the Read and Green's analysis for p-wave superconductors in 
ref.Hl 

can be carried over directly to this model. However, we 
will show that this is not the case and that there are a number 
of subtle differences, the primary one being that our fermions 
do not preserve the gauge symmetries, i.e. the creation of a 
fermion changes the vorticity. 

To proceed we first re-write the fermonic Hamiltonian as 



H 



in 

qq' 



4 °1 



£,qq' &qq' 

a* , -e T , 

qq ^qq 





Cq> 












V 



(28) 



where for example with Hq we would have 

£,qq' 



*2J z&q.q' ~i~ JxX. q(^q,q' —t 



Sq-n :c ,q' 



JyY q(Sq t q' - ny + — n v ,q') 
JxX q (Sq q> — 7 
' q (&q,q f —n 



fiq-n y ,q')- 



(29) 



and the non-zero entries of Hi are given in FIG. [3] To specify 
the particular vortex configuration one, as before, replaces the 
operators X q and Y q by their eigenvalues in that configura- 
tion. On a torus, the £ and A given in d29l ) are modified to 
include the terms that connect both sides of the torus, i.e. the 
terms that connect the sites (0, q y ) to (N x — 1, q y ) and (q x , 0) 
to (q x ,N y — 1). The values of X q and Y q in these terms 
depend on the arrangement of vortices and the quantum num- 



bers 1^ and Iq" of the two independent homologically non- 



,(y) 



(x) (v) 

trivial loop symmetries L and L Q that run through the ori- 
gin (0, 0). The Hamiltonian for any sector on a torus can be 
now be generated by observing the following dependencies 
(see also FIG. g): 



Xq 
Xq 
Xq 

x, 



,q v 
>«» 
,q y 
qy,q* 



11,' =0 



= 1 



— l l In' 



= -I, 



0*0 



w q*,q' y (Qy ¥= and q x ^ N x - 1) 

(q y = and q x ^ N x - 1) 
Wg»,9s (ly ¥= and Qx = N x - 1) 
(q y = and q x = N x - 1) 
(30) 



(X) TT9B-1 



q*,q y 



(.«) 



(q y ±N V -1) 
(Qy = N V -1) 



(31) 



where I 



(y) _ r(v) 



W> nS=o W q ,, qy . These values for X q 

and Y q can be used for the extended Hamiltonian Hi shown 
pictorially in FIG. [3] except for P(q)( 5 > and P(q)( 6 ) when 
q y = N y - 1. For P(q) (5) , we have 



X, 



i = -j£? +I (q*¥=n x -i) 
1), 



Xq x ,q v +1 = Iq X Iq V (qx = N, 



(32) 



while for P(q) (6) 
f X, 



qx,qy 



X, 



q y -l 



9x,9y 



l l 



'W. 



(q x ^N x -1) 
(q x =N X - 1). 



(33) 



The system is diagonalized by solving the Bogoliubov-De 
Gennes eigenvalue problem 





A 




U 


v* 




E 




U 


v* 


At 






V 


u* 




—E 




V 


u* 



t 



(34) 



where the non-zero entries of the diagonal matrix E nm = 
E n Snm are the the quasi-particle excitation energies. The 
Bogoliubov-Valentin quasi-particle excitations are 



ci, cm 



U V* 

V u* 



(35) 
(36) 



which after inversion and substitution into 



give 



M 



1, 



H=Y,Enhh n --). 

71=1 



(37) 



A short review of the relevant theory of the eigenstates and 
eigenvalues of Hamiltonians obtained in this way, in particu- 
lar the application of the Bloch-Messiah-Zumino theorem, is 
given in the appendix. Using the prescriptions described there 
we can calculate the eigenstates and energies in all vortex and 
homology sectors on tori containing a few thousand spins on 
a desktop computer. In the vortex free sector however, there 
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FIG. 4: On a torus there are N/2 + 1 independent loop symmetries^. From these independent loop symmetries all other loop symmetries, 
of which there are 2 JV//2+1 , can be generated by multiplication. We can specify a particular sector of the Hamiltonian by specifying the 
eigenvalues of the N/2 — 1 plaquette symmetries and two homologically non-trivial loops of our choosing. The eigenvalues X q and Y q used 
to specify a particular sector in d28l > are now fully determined. In the graphs above we have indicated how a few specific values of X q and Y q 
depend on the eigenvalues of the independent loop symmetries: (a) Y( 0tNyl ) — —Iq 1 ^, X( Afa ,„ 1 ,o) = —la 3 ) -^(2,3) = I"!, =o W(2,<j y ) ( c ) 

* ( ;w, 2 ) = -#° nj y=0 WW*,*) (d) Y (2iNy _ = 14=0 n^Co 1 w {qmtqv) 



is considerable advantage to be gained by working in the mo- 
mentum representation. The Hamiltonian in this case can be 
written as 

H = £ E k {ili k - (38) 

where the dispersion relation Ej~ is as given for the plane in 
(l24l . The allowed values of k a in the various homology sec- 
tors on the torus are 8 a + 2tt jf- for integer n a =0,1,... N a — 

1, where the four topological sectors, (Iq , Iq ) = (±1, ±1) 

have values of 9 a given by f? Q = {-^-s — )jt- While it is sim- 
ple to see that this expression is valid for the fully periodic 

sector with (1^,1^) = (— 1, —1), it is not so obvious how 
the general expression can be arrived at using the values for 
for X q and Y q given above. One way to understand this jump 
is to imagine that the torus in question is doubled in period 
along the a direction and the pattern of X's and or Y's is re- 
peated on the new part of the lattice, but with a —1 phase. On 
this new lattice it is possible to define a fundamental domain 
of size N x x N y , such that all the values of X q and Y q are 
+1. The periodic boundary conditions on the doubled torus 
now correspond to anti-periodic boundary conditions on the 
fundamental domain and we see that we can use the same dis- 
persion relation as before, only with shifted momenta. 

One would naively expect that the ground-states of the 
Hamiltonian be given by the four analogues of the planar 
ground-state ((27) corresponding to the four ground-states of 
the toric code, and with energies — ^X^fc^fc- These ener- 
gies are not exactly equal for the four homology sectors on 
the torus, because the allowed momenta are different, but it 
is not difficult to see that, at least for the dispersion relations 
we have given, they approach each other rapidly as the system 
size is increased. 

However, there are two situations where the true ground- 
state of a topological and vorticity sector on the torus may not 
be given by the BCS product. First of all, there is a connection 
between the number of vortices and the number of fermions 
on the torus. A configuration with an odd number of vortices 
of electric type can only exist if there are an odd number of 



broken dimers, i.e. an odd number of fermions. BCS products 
like eq. (|27| | have even fermion number parity and hence do 
not apply to vorticity sectors which have an odd number of 
electric (or in fact magnetic) vortices. Secondly, there are sit- 
uations where the allowed momenta together with the values 
of Uk and Vk make the expression for the BCS product state 
vanish. 

A particularly important example of this occurs in the vor- 
tex free sector with (4,4) = (— 1, — 1). In the Abelian 
phase of the model, we can just use the BCS ground-state d27l ) 
as expected. However in the B-phase, we see that A w-ir = 
and ^ 7 r, 7 r/£' 7riir = —1 implying that = and = i. 
This cause the BCS state dZTT i to vanish because on a torus 
c\ v c_ w _ n — (cj. if) 2 — 0. It is important to note that this 
effect is not dependent on the B-phase being gapless, one only 
requires that is negative and A ffi7r = 0. 

The vanishing of the BCS state is somewhat similar to what 
happens in Read and Green's treatment of the spinless p-wave 
superconductor 26 ' 36 , but here, we cannot propose the expres- 
sion in Eq. dA12b as an alternative because it has the wrong 
fermion number parity for the zero vortex sector. However, 
one sees that 

ivo= n (^+v k c{ci k )\{w q }^\^\ {%})., 

(39) 

with {Iq X \ l^) = (—1,-1) and Q q = 1 for all q, is an eigen- 
state with even fermion number, no vortex excitations and en- 
ergy — Ek/2 + -Ett.tt- As with the generic situation shown 
in the Appendix, all of the states | tp^') — 1-n 7r7fc'l VVtt) are 
also vortex-free eigenstates with energy — J^k + Ey- 
The ground-state of the vortex free system in this topologi- 
cal sector is precisely the state in this family for which Ek> is 
minimal. If the system is gapped, then Ek' does not approach 
zero even in the thermodynamic limit, and the ground-state 

of the (?q >4 ) = (~ 1' ~ 1) * s S a PP e d away from the de- 
generate ground-states of the other three vortex free sectors. 
Even more generally, we can say that if the B-phase is gapped 
and the conditions ^^/E^^ = — 1 and A ff . ff = are ful- 
filled, the ground-state on a torus is three-fold degenerate, as 
expected from TQFT. 
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V. SUMMARY AND OUTLOOK 

We have described a spin fermionization procedure for Ki- 
taev's honeycomb model and related spin models. Using this 
method we derived exact expressions for the ground-states 
and their associated eigenvalues. The derived ground-states 
are closed form expressions that do not require additional 
spectral projection. These expression combine two powerful 
wavefunction descriptors: the BCS product and the stabilizer 
formalism. The solution clarifies the nature of the topological 
phases of the model and the role the BCS product plays in de- 
termining them. It is clear now for example that the Z 2 x Z 2 
Abelian phase is determined from the fermionic vacuum, with 
the BCS product only adjusting this state slightly. As the vac- 
uum is fixed and Abelian, the transition to the non-Abelian 
phase is therefore driven exclusively by the BCS product. 

We also showed how to extend our fermionization proce- 
dure to handle general vortex configurations on a torus and 
we discussed how the additional constraints due to the in- 
terdependence of loop symmetries and fermions arise in the 
calculations. We closely examined the ground-state for the 
fully periodic vortex free sector, building on a more general 
discussion for arbitrary configurations given in the Appendix, 
and explained why the blocking mechanisms in the system 
dictate that the non-Abelian ground-state is three-fold degen- 
erate, confirming the prediction from TQFT We intend to use 
our fermionization to study degeneracies in systems with mul- 
tiple anyonic excitations in the future and hope to elucidate the 
relation between the Abelian and non-Abelian anyons in the 
model. 

While we have explicitly derived the relationship between 
the exact solution and the toric code ground-state we have not 
yet explored the relationship between the exact solution and 
its perturbative approximation. This now possible, at least in 
principle, as the BCS product and the Brillouin-Wigner per- 
turbation expansion both start from the TC ground-state. At 
the very least, the comparison should reveal the precise de- 
nominators in the perturbative expansion and clarify how this 
expansion breaks down at the phase transition. More spec- 
ulatively, it may also help to extend perturbative techniques 
beyond the phase transition. The ability to do so could pro- 
vide a new perspective on the non-Abelian phase and would 
be extremely useful for other models which are not exactly 
solvable. 

The results obtained on the torus highlight the connections 
between the blocking mechanisms, the ground-state degener- 
acy and the Abelian to non-Abelian phase transition. A more 
general analysis along these lines, which is not confined to a 
particular model would almost certainly be beneficial. 
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APPENDIX A: GROUND-STATE CONSTRUCTION 

Here we review some of the relevant Hartree-Fock- 
Bogoliubov theory and discuss how to bring the ground-states 
of each vortex sector into a canonical form. We restrict our- 
selves, as in the main text, to situations where we have an even 
number of sites M on the lattice, and therefore an even num- 
ber of 7^ excitations. For more details we refer the reader to 
ref. |43| on which much of the following is based. 

The ground-state for a fermionic Hamiltonian ( f3Tb can usu- 
ally be written down as a Hartree-Fock-Bogoliubov (HFB) 
projection, 

igs)=n^»i vac > (a v 

n 

where the energy of this ground-state is — E n /2. How- 
ever, sometimes the physical situation will demand that the 
ground-state has odd fermion number-parity p. In the honey- 
comb model for example the fermionic number parity is com- 
pletely determined form the vortex configuration. Specifically 
we see that, because creating an e-m vortex pair excitation 
necessary means breaking a z-dimer, the fermionic number 
parity must be equal to the e-number parity and the m-number 
parity. 

The procedure in this case is to redefine our choice of 
7^ and 7 such that the lowest energy odd parity state can 
be found. In practice we swap 7J «-> 71 and then set 
I gs)odd = 7i I gs). This new state, annihilated by the annihila- 
tion operators, has odd fermion number parity but has energy 
— J2 n E n /2 + Ei (Here our convention is to choose E\ to be 
the smallest of the Ei). 

While the above prescription handles a great many physi- 
cal situations there are a number of reasons why it is not al- 
ways sufficient. The first is simply that the method by which 
we construct the ground-state is projective and therefore in 
many cases the vacuum state is not uniquely defined. The sec- 
ond, and perhaps more important, is that physical situations 
do exist such that the calculated 7 fermions are such that | gs) 
and therefore \gs) dd are already zero. A solution to both 
problems can be found by making use of the Bloch-Messiah- 
Zumino theorem2i2£. 

In practice the theorem says that we can do a singular value 
decomposition of the M x M matrix U 

U = DUC (A2) 

for unitary C and D such that the eigenvector matrix defining 
the 7T quasi-particle excitations, d35l l. can be decomposed into 



u 


v* 




D 







u 


V 




c 





V 


u* 







D* 




V 


u 







c* 



(A3) 



where U and V take block diagonal forms 
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The ground-state, written in the canonical basis, is 



U = 



(A4) 



and 



V 



Vi 



(A5) 



where Z and / are square zero and identity matrices respec- 
tively and 



Ui 



and 



m 
Ui 



v % 
-vt 



(A6) 



(A7) 



This factorization means that d35l l can be understood as 
three separate transformations: 

(1) A unitary operator that mixes the fermionic excitation and 
annihilation operators amongst themselves: 



(A8) 



3 J 

This defines what is known as the canonical basis. 



(2) A Bogoliubov Transform which, for paired levels (u k > 
0, v k > 0), mixes creation and annihilation operators 

a k = u k a\ - v k a- k a k = u k a\ + v k a k (A9) 



where the (k, k) label elements of the 2x2 matrices 
and dA7b . and for blocked levels which are either occupied 

(Vi = 1, u-i = 0) or empty (vj = 0, Uj = 1): 



at = a,i, a] = a], 
Ui — aj, ctj = cij. 



(A10) 



(3) A unitary operator which mixes quasi-particle operators 
a\. amongst themselves: 



7/ = Ckia l ^ 



k 



(All) 



I s ) 



ihiik 

i k^i 



v k a\a 



|)|vac), 



(A12) 



where the first product is over the m occupied levels only and 
it is understood that all m's in the second product are non- 
zero. The state is annihilated by all j n and thus has an energy 
— J2 n En/2- To see this first note that the state is annihilated 
by all cti and then that each 7; is a linear superposition of 
these. 

The fermion number-parity p of the state is dictated by the 
number of occupied modes m. The appearance of an odd 
number of occupied modes implies that the final state is of 
different number parity to the vacuum. This is exactly the situ- 
ation observed by Read and Green for the spinless p-wave su- 
perconductor on a toru s - 26 ' 36 . In the honeycomb lattice model 
on a torus, the situation is complicated slightly by the fact 
that the fermion number parity of an eigenstate is determined 
exclusively by the vortex configuration sector to which it be- 
longs. This immediately implies that an odd (even) number 
of occupied modes is not allowed in a vortex sector with even 
(odd) e-number or m-number parity. When these situations 
occur we say that the state (IA12b is blocked and we must 
re-arrange the eigenvector matrix, effectively switching the 
7 J and 7], such that the eigenstate with the lowest energy is 
achieved. Below we discuss t he four possible scenarios and 
describe how to construct the eigenstate in each case. 

(A) Even p, Even m: 

Here we may use ( IA12I ) with no modification. The ground- 
state energy is E m [ n = — E n /2. One encounters this sit- 
uation in the vortex free sector in both the A-phase and the 
three partially/fully anti-periodic sectors of the B-phase. 

(B) Even p, Odd m: 



Here we also use (IA12I) but the singular value decompo- 
sition to calculate the u's and w's is performed after first 
switching columns (Un,Vn) <-> (VJ*,?/^). Of course an 
initial singular value decomposition of the original U matrix 
is first needed in order to determine the number of occupied 
modes m. The switching of columns of the matrix effectively 
changes an occupied mode for an empty one and the energy of 
the ground-state is therefore E m i n — — J2 n E n /2 + E%. We 
see that if the sector has even fermion parity and is gapped, 
an odd number of occupied modes rn implies a raising of 
the energy above what one might otherwise expect. We en- 
counter this situation in the B-phase of the vortex-free sector 
with (/<»>,#>) = (-1,-1). ' 

(C) Odd p, Even m: 

This situation is again handled by switching occupied and 
empty modes (Un, Vn) <-> (VJ*, U n ) and the energy of the 
ground-state is again E m [ n = — J2 n E n /2 + E\. The vacuum 
in this case must be from the vortex sector such that the op- 
eration of the now odd number of aj's gives the vortex and 
topological sector for which we calculated the U and V ma- 
trices. 
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(D) Odd p, Odd to. 

The ground-state in this sector is given by (IA12I) with- 
out modification and the ground-state energy is — J2 n E n /2. 
Thus, if the sector is gapped, the energy of this ground-state is 
lower than what one might expect from an odd fermion num- 
ber state. As in the previous case the vacuum must be defined 
so that operating with an odd number of aj's gives the correct 
vortex and topological sector. 

In all of the above situations, because the jj's do not com- 



mute with all loop symmetries, the excited states in each 
sector are obtained by operating on the ground-state with 
quadratic operators 'yjjj. This can also be checked through 
a simple counting argument. On a torus, as we have M + 1 
independent loop symmetries, the Hilbert space dimension of 
each sector is 2 M_1 , see- . If we were allowed to operate with 
single 7^'s and 7/s on the ground-state we could generate 2 M 
states, which is obviously too many. 
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